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ABSTRACT 

The anisotropy of clustering in redshift space provides a direct measure of the growth 
rate of large scale structure in the Universe. Future galaxy redshift surveys will make 
^_ ) . high precision measurements of these distortions, and will potentially allow us to dis- 

tinguish between different scenarios for the accelerating expansion of the Universe. 
Accurate predictions are needed in order to distinguish between competing cosmolog- 
ical models. We study the distortions in the redshift space power spectrum in ACDM 
and quintessence dark energy models, using large volume N-body simulations, and 
test predictions for the form of the redshift space distortions. We find that the lin- 
ear perturbation theory prediction by Kaiser (1987) is a poor fit to the measured 
distortions, even on surprisingly large scales k fJ.05ft.Mpc~ 1 . An improved model 
for the redshift space power spectrum, including the non-linear velocity divergence 
CN ' power spectrum, is presented and agrees with the power spectra measured from the 

simulations up to k ~ 0.2/iMpc -1 . We have found a density-velocity relation which is 
cosmology independent and which relates the non-linear velocity divergence spectrum 
to the non-linear matter power spectrum. We provide a formula which generates the 
non-linear velocity divergence P(k) at any redshift, using only the non-linear matter 
power spectrum and the linear growth factor at the desired redshift. This formula is 
CO accurate to better than 5% on scales k < 0.2/iMpc -1 for all the cosmological models 

discussed in this paper. Our results will extend the statistical power of future galaxy 
surveys. 
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1 INTRODUCTION as ESA's Euclid mission (jCimatti et alj|2009h . will be able 

to achieve this precision. These redshift space distortions 
The rate at which cosmic structures grow is set by a competi- , , „ , ,. . , .. ,, 

are commonly modelled using a linear perturbation theory 
tion between gravitational instability and the rate of expan- . TT , . . ., ,. ... „,,. 

, . , expression. We test the validity ot this approximation using 

sion of the Universe. The growth of structure can be mea- , , , T , , . , , . , , , , , . 

, .. large volume N-body simulations to model the redshift space 

sured by analysing the distortions m the galaxy clustering ,. . ,. . »^,t^iv«- i ■ , , i i i 

distortions m ACDM and quintessence dark energy models, 
pattern, when viewed m redshift space i.e. when a galaxy s ., ,,.111 ■ j , , , j r 

v to see it it works at the level required to take advantage of 

redshift is used to infer its radial position). Proof of con- . „ ,. . r . „, , , „ 

f r-— - ' — : rn 1 -j the information m forthcoming surveys. 1 he large volume of 

cept of this approach came recently from Guzzo ct al. ( 2008 . . .. . . „ , ,. ., r 

: ~ ' , ' our simulations means that we are able to hnd the limits of 

who used spectroscopic data for 10,000 galaxies from the , .. , , ... , , 

. t ~n t 1 — 1 . . t I ^ — ^ ' ~i 1 1, perturbation theory models. We can also study the impact 

VIMOS-VLT Deep Survey the Fevre et all 12005) to mea- F , .. ... y , , . . . ... * 

, ~ — r " of non-lmeanties on large scales m cosmologies with differ- 

sure the growth rate of structure at redshift z = 0.77 to an ,. , . . j™,, , . . 

rr : -11 1 cut expansion histories from ACDM, such as quintessence 

accuracy of ~ 40% (see also Peacock ct al. 2001). To distin- , , 

— J r 1 — dark energy. 

guish between competing explanations for the accelerating 

expansion of the Universe, we need to measure the growth of One explanation of the accelerating expansion of the 

structure to an accuracy of a few percent over a wide redshift Universe is that a negative pressure dark energy component 
interval. The next generation of galaxy redshift surveys, such makes up approximately 70% of the present densit y of the 

Universe (|Komatsu et a l. 2009; Sanchez ct al. 200^). Exam- 
ples of dark energy models include the cosmological constant 
* E-mail: elise.jennings@durham.ac.uk and a dynamical scalar field such as quintessence ( see e.g. 
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ICopeland et~aHl2006L for a review). Other possible solutions 
require modifications to general relativity and include exten- 
sions to the Einstein-Hilbert acti on, such as f(R) theories 
or br aneworld cosmologies (see e.g. lDvali et al.ll2000l ; IOvaizul 
l200Sh . 

The expansion history of the Universe is described by 
the scale factor, a(t). Dark energy and modified gravity mod- 
els can produce similar expansion histories for the Universe, 
which can be derived from the Hubble parameter measured, 
for example, using type la supernovae. As both dark energy 
and modified gravity models can be described using an effec- 
tive equation of state which specifies the expansion history, 
it is not possible to distinguish between these two possibili- 
ties using measurements of the expansion history alone. 

The growth rate is a measure of how rapidly struc- 
tures are forming in the Universe. Dark energy or modified 
gravity models predict different growth rates for the large 
scale structure of the Universe, which can be measure d using 
redshi ft space distortions of clustering. As noted by ILinderl 
(2005), in the case of general relativity, the second order 
differential equation for the growth of density perturbations 
depends only on the expansion history through the Hub- 
ble parameter, H(a), or the equation of state, w(a). This is 
not the case for modified gravity theories. By comparing the 
cosmic expansion history with the growth of structure, it is 
possible to distinguish the physical origin of the accelerat- 
ing expansion of the Univer se as being due either to dark 
energy or modified gravity l|Lue. Scoccimarro fc Starkmanl 
12004 : ; ILinderl 120051 ). If there is no discrepancy between the 
observed growth rate and the theoretical prediction assum- 
ing general relativity, this implies that a dark energy com- 
ponent alone can explain the accelerated expansion. 

Galaxy redshift surveys allow us to study the 3D spatial 
distribution of galaxies and clusters. In a homogeneous uni- 
verse, redshift measurements would probe only the Hubble 
flow and would provide accurate radial distances for galax- 
ies. In reality, peculiar velocities are gravitationally induced 
by inh omogeneous st ructure and distort the measured dis- 
tances. iKaiserl (|l987T ) described the anisotropy of the clus- 
tering pattern in redshift space but restricted his calculation 
to large scales where linear perturbation theory should be 
applicable. In the linear regime, the matter power spectrum 
in redshift space is a function of the power spectrum in real 
space and the parameter f3 — f /b where / is the linear 
growth rate. The linear bias factor, b, characterises the clus- 
tering of gal axies with re s pect to the underlyi ng mass distri- 
bution (e.g. Kaiser 1987). S coccimarr o (2004) extended the 
analysis of Kaiser (1987) into the non-linear regime, includ- 
ing the contribution of peculiar velocities on small scales. 
We test this model in this paper. 

Perturbations in bulk flows converge more slowly then 
perturbations in density, and so very large volume simula- 
tions are needed to model these flows, and hence the redshift 
space distortion of clustering, accurately. Our simulation 
boxes are 125 times the volume of those used bv lCole et al.l 
1 19941 ) and ~ 30 times the volum e of the N-body res u lts in - 
terpreted by IScoccimarrol (12004) 1. IPercival fc White! (|2009l ) 
used a single l/i _1 Gpc box to study redshift space distor- 
tions in a ACDM model. Their simulation is over three time 
smaller than the one we consider. 

This paper is organised as follows: In Section [2] we dis- 
cuss the linear growth rate and review the theory of redshift 



space distortions on linear and non-linear scales. In Section 
[3] we present the quintessence models considered and the 
details of our N-body simulations. The main results of the 
paper are presented in Sections [4] and [5] The linear the- 
ory redshift space distortion, as well as models for the red- 
shift space power spectrum which include non-linear effects 
are examined in Section [4] for various dark energy cosmolo- 
gies. In Section [5] we present the density- velocity relation 
measured from the simulations. Using this relation the non- 
linear models used in the previous section can be made cos- 
mology independent. We present a prescription for obtaining 
the non-linear velocity divergence power spectrum from the 
non-linear matter power spectrum at an arbitrary redshift 
in Section \5. 21 Our conclusions are presented in Section [6] 



2 REDSHIFT SPACE DISTORTIONS 

In Section 12.11 we consider several parametrizations which 
are commonly used for the linear growth rate. In Section 
12.21 we review linear perturbation theory for redshift space 
distortions and discuss the assumptions that are used in this 
approach. In Section ^. 3l we present several models proposed 
to describe the distortio ns in the non-linea r regim e. A similar 
review can be found in IPercival fc White! ([20091 ). 



2.1 Linear growth rate as a probe of gravity 



The linear gro wth rate is a promising probe of the nature 



of dark energy jGuzzo et al.l 2008: Wang 


120081: Linderl 20081: 


Song 


fc Percival 


200£ 


: White et al.l 2009 


: Percival & White! 


2009 


; IStril et al. 


|2010 


;ISimDSon & PeacockH2010l). Although 



the growth equation for dark matter perturbations is easy to 
solve exactly, it is common to consider parametrizations for 
the linear growth rate, / = dlnD/dlna, where D(a) is the 
linear growth factor. These parametrizations employ differ- 
ent variables with distinct dependencies on the expansion 
and growth histories. 

A widely used approximatio n for /, first sugge sted by 
iPeebles! l| 19761 ), is f(z) « fl^ 6 . lLahav et all (|l99lh found 
an expression for /, in terms of the present day densi- 
ties of matter, S7 m , and dark energy, f2 DE , which showed 
only a weak dependence on th e dark e nergy density, with 
/ « Q.™ + n DE /70(l + n m /2) ILinderl (120051 ) extended the 
analysis of I Wang fc Steinhardtl 1 19981 ) to find a new fitting 
formula to the exact solution for the growth factor, which 
he cast in the following form 

g(a) = ^ a exp Qf* dlna [nj(o) - 1]) , (1) 

where 7 is the index which parametrises the growth his- 
tory, while the expansion his tory is described by the matter 
density Q m (a). ILinderl (|2005l ) proposed the empirical result 
7 = 0.55 + 0.05[1 + w(z = 1)], where w is the dark energy 
equation of state, w hich gives / = tljj; 55 for a cosmologi- 
cal constant (see also lLinder fc Cahnll2007l ) . We discuss this 
formula for / further in Section [3] when we introduce the 
quintessence dark energy models used in this paper. 



Redshift space distortions 3 



2.2 Linear redshift space distortions 

The comoving distance to a galaxy, s, differs from its true 
distance, x, due to its peculiar velocity, v(x) (i.e. an addi- 
tional velocity to the Hubble flow), as 



x + 



H (a) 



(2) 



where H(a) is the Hubble parameter and v ■ x, is the pecu- 
liar velocity along the line of sight. Inhomogeneous structure 
in the universe induces peculiar motions which distorts the 
clustering pattern measured in redshift space on all scales. 
This effect must be taken into account when analyzing three 
dimensional datasets which use redshift as the radial coor- 
dinate. Redshift space effects alter the appearance of the 
clustering of matter, and together with non-linear evolution 
and bias, lead the power spectrum to depart from simple 
linear perturbation theory predictions. 

On small scales, randomised velocities associated with 
viralised structures decrease the power. The dense central re- 
gions of galaxy clusters look elongated along the lin e of sight 
in re dshift space, which produces 'fingers of God' (|Jacksonl 
1 19721 ) in redshift survey cone plots. On large scale s, coherent 
bulk flows distort clustering statistics, (see H amilton! [l998l . 
for a review of redshift space distortions). For growing per- 
turbations on large scales, the overall effect of redshift space 
distortions is to enhance the clustering amplitude. Any dif- 
ference in the velocity field due to mass flowing from under- 
dense regions to high density regions will alter the volume 
element, causing an enhancement of the apparent density 
contrast in redshift space, 8 3 (f), compared t o that i n rea l 
space, & r (r). This effect was first analyzed bv lKaiseil (|l987f ) 
and can be approximated by 



(3) 



where fi is the cosine of the angle between the wavevector, 
fc, and the line of sight, f3 = f /b and the bias, 6=1 for dark 
matter. 

The Kaiser formula (Eq. [3]) relates the overdensity in 
redshift space to the corresponding value in real space using 
several approximations: 

1. The small scale velocity dispersion can be neglected. 

2. The velocity gradient |du/dr| <€. 1. 

3. The velocity and density perturbations satisfy the 
linear continuity equation. 

4. The real space density perturbation is assumed to 
be small, \S(r)\ <C 1, so that higher order terms can be 
neglected. 

All of these assumptions are valid on scales that are well 
within the linear regime and will break down on different 
scales as the density fluctuations grow. The linear regime 
is therefore defined over a different range of scales for each 
effect. 

The matter power spectrum in redshift space can be 
decomposed into multipole moments using Legendre poly- 
nomials, Li(fi), 



(4) 



The anisotropy in P(k) is symmetric in fj,, as P(k, fi) = 
P(k, — n), so only even values of I are summed over. Each 



multipole moment is given by 

91 -+- 1 r 1 
Pi{k) = ^— / P(fe )M )L i (M)dM 1 



(5) 



where the first two non-zero moments have Legendre poly- 
nomials, Lo(fi) — 1 and I/2(/i) = (3/i 2 — l)/2. Using the 
redshift space density contrast, Eq. [3] can be used to form 
P(k, /i) and then integrating over the cosine of the angle /i 
gives the spherically averaged monopole power spectrum in 
redshift space, Pg(k), 



(6) 



P''(fc) 3 7 5 7 ' 

where P r (k) denotes the matter power spectrum in real 
space. In practice, P r (k) cannot be obtained directly 
for a real survey witho ut making approximations (e.g. 
iBaugh fc EfstathioJ[l994T ). 

In th is paper we also consider the estimator for / sug- 
gested bv lCole et ail l|l994h . which is the ratio of quadrupole 
to monopole moments of the redshift space power spectrum, 
Pi(k)/Po(k). From Eq. [3] and after spherically averaging, 
the estimator for / is then 



Pi(k) _ 4//3 + 4/ 2 /7 



Pi(k) 1 + 2//3 + /V5 



(7) 



which is independent of the real space power spectrum. Here, 
as before, / = j3/b, with 6=1 for dark matter. 



2.3 Modelling non-linear distortions to the power 
spectrum in redshift space 

Assuming the line of sight component is along the z-axis, 
the fully non-linear relation between th e real and redshift 
space power spectrum can be written as (|Scoccimarro et all 
1999) 



S(fc '" ) = /(!r)3 



r <e lAA "^(x)-/V z -Mx)] 

x[5(x')-/V;-«*(x')]>, (8) 



where A = fk/i, u z is the comoving peculiar velocity along 
the line of sight, Au z = u z (x) — u z (x'), r — x — x' and the 
only approximation made is the plane parallel approxima- 
tion. This expression is th e Fourier analog of the 'streaming 
model' first sugg ested by iPeebles! jl98d) and modified by 
iFishej (jl995h to take into account the density-velocity cou- 
pling. At small scales (as k increases) the exponential com- 
ponent damps the power, representing the impact of ran- 
domised velocities inside gravitationally bound structures. 

Simplified models for redshift space distortions are 
frequently used. Examples include multiplying Eq. [6] by 
a factor which attempts to take into account small 
scale effects and is eith er a Gaussian or an exponential 
ijPeacock fe Doddsl ll994). A popular phenomenological ex- 
ample of this which incorporates the damping effect of ve- 
locity di spersion on small scales is the so called 'dispersion 
model' (|Peacock fc Dodds!ll994l '). 



P s (k,n)=P r (k)(l + (3n 2 f 



(l + fcV<^/2) 



(9) 



where a v is the pairwise velocity dispersion along the line of 
sight, which is treated as a parameter to be fitted to the 
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Figure 1. Left panel: The linear growth factor divided by the scale factor as a function of redshift for the SUGRA and CNR quintessence 
models and ACDM, as indicated by the key. Right panel: The linear growth rate, / = dlnD/dlna, for the two dark energy models and 
ACDM as a function of redshift. In both the left and right main panels, solid lines represent the exact solution for the linear growth 
factor and growth rate and dashed lines show the fitting formula given in Eq. [T] Note in the right main panel th e ACDM grey dashed 
line has been omitted for clarity. The lower left hand panel shows the formula for D(a)/a given by Lindcr (2005) divided by the exact 
solution as a function of redshift. The ratio of the formula in Eq. [T] for the growth rate, /, to the exact solution is shown in the lower 
right hand panel. Also in the lower right panel the dotted lines show the ratio of the fitting formula / = fl^ 6 to the exact solution for 
each of the dark energy models plotted as a function of redshift. 



data. Using numerical simulations, lHatton fc Cole] <|l999h 
found a fit to the quadrupole to monopole ratio Pi/Po = 
(P£/Po )iin(l — x 1 ' 22 ) to mimic damping and non-linear ef- 
fects, where [P% /Po)ua is the linear theory prediction given 
by Eq.0 x = k/k\ and fci is a free parameter. They extended 
the dynamic range of simulations, to replicate the effect of 
a larger box, using the approxi mate method for adding long 
wavelength power suggested bv lColel (|!997l 1. 

The velocity divergence auto power spectrum is the em- 
semble average, Pgg = (|#| 2 ) where 8 = V • u is the velocity 
divergence. The cross power spectrum of the velocity di- 
vergence and matter density is Pse — (\S8\), where in this 



notation the matter density auto spectrum is Ps, 



2 >- 



In Eq. [8] the term in square brackets can be re-written in 
terms of these non-linear velocity divergence power spec- 
tra by multiply ing out the b r acket s and using the fact that 
Hi = ki ■ z/ki. IScoccimarrol (2004) proposed the following 
model for the redshift space power spectrum in terms of 
Pss, the non-linear matter power spectrum, Pee and Pse, 

P s (fc,/x) = (10) 
(Pss(k) + 2ffi 2 Pse(k) + ffi 4 Pee(k)) x e~ (Sk ^? , 

where a v is the ID linear velocity dispersion given by 



2 1 f Peejk) j3u 



(11) 



In linear theory, Pee and Pse take the same form as Pss and 
depart from this at different scales. Using a simulation with 



512 3 particles in a box of length 479/i _1 Mpc (jYoshida et all 
1200 if ) . IScoccima rro (2004) showed that this simple ansatz for 
P s (fc,/x) was an improvement over the Kaiser formula when 
comparing to N-body simulations in a ACDM cosmology. As 
this is a much smaller simulation volume than the one we 
use to investigate redshift space distortions we are able to 
test the fit to the measured power spectrum on much larger 
scales and to higher accuracy. 



3 N-BODY SIMULATIONS OF DARK 
ENERGY 

In the following sections we briefly review the quintessence 
models discussed in this paper and the N-body simulations 
used to measure various power spectra. 



3.1 Quintessence models 

In quintessence models of dark energy, the cosmological 
constant is replaced b y an extremely light s calar field 
which evolves slowly jRatra fc Peebles! 1 19881; | Wetterichl 



1 19881 ; ICaldwell et aT1ll99Sl ; iFerreira fc Jovcel ll998T ). Differ^ 
ent quintessence dark energy models have different dark 
energy densities as a function of time, Q-d E (z). This im- 
plies a different growth history for dark matter perturba- 
tions from that expected in ACDM. In this paper we con- 
sider three quintessence models, each with a different evolu- 
tion for the dark energy equation of state parameter, w(a). 
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Figure 2. Left panel: The ratio of the monopole redshift power spectra and real space power spectra measured from the ACDM simulation 
at z = and z = 1 are plotted as blue lines. The error bars plotted represent the scatter between the different power spectra from four 
ACDM simulations set up with different realisations of the density field with the distortions imposed along either the i,j or z axis and 
averaged. The power spectra P(fc,/i = kx/k), P(k,/i = k y /k) and P(k,/i = k z /k) measured from one simulation are plotted as the cyan, 
purple and red dashed lines respectively. Right panel: The ratio of the quadrupole to monopole moment of the redshift space power 
spectrum measured from the simulations at z = and z = 1 in ACDM are plotted in blue. It was not possible to accurately measure 
the quadrupole to monopole power in the first bin, so this point has not been plotted in the right hand panel. Note for wavenumbers 
k > 0.1/iMpc -1 , only every fifth error bar is plotted for clarity. The Kaiser formula, given by Eq. [B] is plotted as a blue dotted line. The 
error bars were obtained as described for the left-hand panel. 



These models are a representative sample of a range of 
quintessence models a nd are a subset of those considered by 
IJennines et"aH (|201fjh to which we refer the reader for fur- 
ther d etails. Briefly, the SUGRA model of iBrax fc Martini 
( 1999) has an equation of state today of Wo = —0.82 and a 
linear growth factor which differs from ACDM by 20% at 
z = 5. The 2EXP model has an equation of state that makes 
a rapid transition to too = —1 at z = 4 and since then has a 
similar expansion history to ACDM ijBarreiro et al.ll2000h . 
The CNR quintessence model has a non- negligible amount 
of dark energy at early times and an eq uation of state to- 
day of too = —1 (|Copeland et alj|2000l 1. The dark energy 
equation of state for each model is described using a 4 vari- 
able parametrization for w(a) which is able to accurately 
describe the expansion history ove r the full range of red- 
shifts modelled by the simulations |Corasaniti k, Copelandl 
120031 ). 

The presence of small but appreciable amounts of dark 
energy at early times also modifies the growth rate of fluc- 
tuations from that expected in a matter dominated uni- 
verse and hence changes th e shape of the linear t heory P(k) 
from the ACDM prediction l| Jennings et al.ll2010f ). The CNR 
quintessence model used in this paper has non-negligible 
amounts of dark energy at high redshift s and so could be 
classe d as an 'early dark energy' model (|Doran fc Robberi 
2006). As a result, the linear theory power spectrum is ap- 
preciably different from that in a ACDM co smology, with 
a broader turnover, (see IJennines et~aT1l2010l , for futher de- 
tails). 

Quintessence dark energy models will not necessarily 
agree with observational data if we adopt the same cosmo- 
logical parameters as used in the best fitting ACDM cos- 
mology. These best fit parameters were found using the ob- 
servational constraints on distances such as the angular di- 



Table 1. Cosmological parameters used in the simulations. The 
first column gives the cosmological model, the second the present 
day matter density, O m , the third the baryon density, Q^, and the 
fourth the Hubble constant, h, in units of 100 km s _1 Mpc . 



Model 


S2m 


n b 


h 


ACDM / 2EXP 


0.26 


0.044 


0.715 


SUGRA 


0.24 


0.058 


0.676 


CNR 


0.28 


0.042 


0.701 



ameter distance to last scattering and the sound horizon at 
this epoch, from the cosmic microwave background, as well 
as distance measurements from t he baryonic acoustic o scil- 
lations and Type la supernovae (|jennings et alj l2010h . In 
this paper the best fitting cosmological parameters for each 
quintessence model are used in the N-body simulations, as 
listed in Table 1. 

In the left panel of Fig. [1] we plot the exact solution for 
the linear theory growth factor, divided by the scale factor, 
as a function of redshift together with the fitting formula in 
Eq. [T] The 2EXP quintessence model is not plotted in Fig. [T] 
as the linear growth factor for this model diff ers from ACDM 
only at high redshifts, z > 10. iLinderl (|2005h found that the 
formula in Eq.[T]reproduces the growth factor to better than 
0.05% for ACDM cosmologies and to ~ 0.25% for different 
dynamical quintessence models to the ones considered in 
this paper. We have verified that this fitting formula for D 
is accurate to ~ 1% for the SUGRA and 2EXP dark energy 
models used in this paper, over a range of redshifts. Note, in 
cosmological models which feature non negligible amounts 
of dark energy at high redshifts, a fu rther correct ion fac- 
tor is needed to this parametrisation (lLinderll2009l), Using 
the parametrization for w(a) provided by Doran fc Robbera 
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(2006) for 'early dark energy ', lLindeij f2009) proposed a sin- 
gle correction factor which was independent of redshift. The 
CNR model has a high fractional dark energy density at 
early times and as a result we do not expect the linear the- 
ory growth to be accurately reproduced by Eq. [1] As can 
be seen in Fig. [T] for the CNR model, any correcti o n fac - 
tor between the fitting formula suggested bv lLinderl (|2005T l 
and the exact solution for D/a would depend on redshift 
and is not simply a const ant. In this case, t h e 'ear ly dark 
energy' parametrisation of iDoran fc RobbersI (|2006h is not 
accurate enough to fully describe the dynamics of the CNR 
quintessence model. This difference is ~5% at z — 8 for the 
CNR model, as can be seen in the ratio plot in the left panel 
of Fig [T] The exact solution for the linear growth rate, /, 
and the fitting formula in Eq. [T] / = fl^ (a), is plotted in 
the right panel of Fig. [I] The old approximation / = ^f', is 
plotted in the bottom right panel in Fig. [T] The dotted lines 
represent the ratio / = fi^ 6 to the exact solution for each of 
the dark energy models. It is clear that this approximation 
for the growth factor is not as accurate as the formula in 
Eq. [T]over the same range of redshifts. 

3.2 Simulation Details 

We use the N-b ody simulations carried out by 
IJennings et"ai] (|201Ch . These simulations were performed at 
the Institute of Computational Cosmology using a memory 
efficient versio n of the Tree PM code Gadget-2, called 
L-Gadget-2 (|Springell [20051 ). For the ACDM model we 
used the following cosmological parameters: f2 m = 0.26, 
SIde = 0.74, n b = 0.044, h = .715 and a spectral tilt 
of ns = 0.96 (jSanchez et al.| [2009). The linear theory rms 
fluctuation in spheres of radius 8 Mpc is set to be 

erg = 0.8. For each of the quintessence models, a four 
variable parametrization of the dark energy equation of 
state is used as described above. In each case, the cos- 
mological parameters used are the best fitting parameters 
to observational constraints from the cosmic microwave 
background, baryonic acoustic oscillations and supernovae 
la taking into account the impact of the quintessence m odel. 
(Stage III in the terminology of I Jennings et alj|2~010h . 

The simulations use N = 646 3 ~ 269 x 10 6 particles 
to represent the matter distribution in a computational box 
of comoving length 1500/i _1 Mpc. The comoving softening 
length is 50h~ kpc. The particle mass in the ACDM simu- 
lation is 9.02 x 10 11 /i _1 Mq and is slightly different in the 
other runs due to changes in fi m (see Table 1). The initial 
condition s were set up starting from a gla ss configuration of 
particles l|Whitelll994l ; iBaugh et alJll995h . In order to limit 
the impact of the initial displacement scheme we chose a 
starting redshift of z = 200. 

The linear theory power spectrum used to gener- 
ate the initial con ditions was obtained using CAMB 
|Lewis fc Bridldl2002h . We use a modified version of CAMB 
which incorporates the influe nce of dark energy on dark mat- 
ter clustering at early times (|Fang et al.ll2008l ). 

In each model the power spectrum at redshift zero is 
normalised to have erg = 0.8. Using the linear growth fac- 
tor for each dark energy model, the linear theory P(k) 
was then evolved backwards to the starting redshift of 
z — 200 in order to generate the initial conditions. The 
power spectrum was computed by assigning the particles 



to a mesh using the cloud in cell (CIC) assignment scheme 
|Hocknev fc Eastwood! [19811 ) and performing a fast Fourier 
transform (FFT) of the density field. To compensate for the 
mass assignment sch eme we perform an ap proximate de- 
convolution following iBaumgart fc Frvl (|l99ll ). 



4 RESULTS I: THE MATTER POWER 
SPECTRUM IN REAL AND REDSHIFT 
SPACE 

In Sections 14.11 and 14.21 we present the redshift space dis- 
tortions measured from the simulations in ACDM and 
quintessence cosmologies, and we compare with the predic- 
tions of the linear and non-linear models discussed in Section 

4.1 Testing the linear theory redshift space 
distortion 

In the left panel of Fig. we plot the ratio of the redshift 
space to real space power spectra, measured from the ACDM 
simulation at z — and z = 1. Using the plane parallel ap- 
proximation, we assume the observer is at infinity and as a 
result the velocity distortions are imposed along one direc- 
tion in fc-space. If we choose the line of sight direction to 
be the z-axis, for example, then /x = k z /k where k = |fc|. 
In this paper the power spectrum in redshift space repre- 
sents the average of P(k,/j, = k x /k), P(k,fi = k y /k) and 
P(k,n = k z /k) where the line of sight components are par- 
allel to the x, y and z directions respectively. We use this 
average as there is a significant scatter in the amplitudes 
of the three redshift space power spectra on large scales, 
even for a computational box as large as the one we have 
used. The three monopoles of the redshift space power spec- 
tra P(k,n = k x /k), P(k,n = ky/k) and P(k,/i = k z /k) 
measured in one of the realisations are plotted as the cyan, 
purple and red dashed lines respectively, to illustrate the 
scatter. 

In Fig.[2]the Kaiser formula, given by Eq.|6] is plotted as 
a blue dotted line, using a value of / = 0^ 55 (z) for ACDM. 
The error bars plotted represent the scatter over four real- 
isations after averaging over P(k) obtained by treating the 
x,y and z directions as the line of sight. It is clear from 
this plot that the linear perturbation theory limit is only at- 
tained on extremely large scales (k < 0.03/iMpc _ ) at z = 
and at z = 1. Non- linear effects are significant on scales 
0.03 < fc(/iMpc _1 ) < 0.1 which are usually considered to be 
in the linear regime. The measured variance in the matter 
power spectrum on these scales is 10~ 3 < a 2 < 10 -2 . 

In the right panel of Fig. 2 we plot the ratio Pi/Po 
for ACDM at z — and 2 = 1. The ratio agrees with 
the Kaiser limit (given in Eq. 7) down to smaller scales, 
k < 0.06/iMpc" 1 , compared to the monopole ratio plotted 
in the left panel. Our results agree with previous work on 
the quadrupole and monop ole moments of the redshift space 
power spectrum for ACD M (|Cole et al.1 19941 ; lHatton fc Cold 
1 19991 ; IScoccimarrd 120041 ) . At z = 1, the damping effects are 
less prominent and the Kaiser limit is attained over a slightly 
wider range of scales, k < 0.1/iMpc _1 , as non-linear effects 
are smaller then at z = 0. In the next section, we consider 
these ratios for the quintessence dark energy models in more 
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Figure 3. Left panel: The ratio of the non-linear power spectra, Pgg, Pse an d Peg for ACDM measured from the simulation at z = 0, 
divided by the corresponding power spectrum measured from the simulation at 2 = 5, scaled using the square of the ratio of the linear 
growth factor at z = 5 and 2 = 0. The non-linear matter power spectrum is plotted as a grey dot-dashed line, the non-linear velocity 
divergence auto power spectrum Pee is plotted as a blue solid line and the non-linear cross power spectrum, Pse, is plotted as a green 
dashed line. Right panel: The ratio of the non-linear power spectra, Pss, Pse an d Pee? to the linear theory matter P(k) in ACDM 
measured from the simulation at z = 0. All power spectra have been divided by the linear theory matter power spectrum measured from 
the simulation at z = 5, scaled using the square of the ratio of the linear growth factor at z = 5 and z = 0. In both panels the error bars 
represent the scatter over eight ACDM realisations after imposing the peculiar velocity distortion along each cartesian axis in turn. 
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Figure 4. Left panel: The ratio of the non-linear power spectra, Pss, Pse an< i Pee? to the linear theory P(k) in ACDM measured from 
one realisation of the matter density and velocity fields at z = 0. All power spectra have been divided by the linear theory matter power 
spectrum measured from the simulation at z = 5, scaled using the square of the ratio of the linear growth factor at z = 5 and z = 0. 
Right panel: Similar to that in the left panel but for the SUGRA quintessence model. The lines are the same as used in the left hand 
panel. 



detail. For each model we find that the analytic expression 
for the quadrupole to monopole ratio describes the simu- 
lation results over a wider range of wavenumber then the 
analogous result for the monopole moment. 



4.2 Nonlinear models of P s (k,fl) 

The linear theory relationship between the real and redshift 
space power spectra given in Eq. \§\ assumes various non- 
linear effects are small and can be neglected on large scales. 
These assumptions are listed in Section \2.2\ In this section 



we consider the non-linear terms in the gradient of the line 
of sight velocity field and explore the scales at which it is 
correct to ignore such effects in the redshift space power 
spectrum. As a first step, we compare the model in Eq. 
1 1 01 to measurements from N-body simulations for different 
quintessence dark energy models, without the damping term 
due to velocity dispersion. This will highlight the scale at 
which non-linear velocity divergence terms affect the mat- 
ter power spectrum in redshift space and cause it to depart 
from the linear theory prediction. 

If we rewrite dS/dr as aH(a)f(Q m (a),j) 5, where S is 
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Figure 5. A comparison of the impact of the FFT grid dimen- 
sion on power spectrum estimation. The plots show the ratio of 
the non-linear power spectra, Pgg (upper panel) and Pgg (lower 
panel), to the linear theory matter power spectrum measured 
from the simulations in ACDM, using different FFT grid sizes. 
From bottom to top in each panel the lines show the ratios for 
grid sizes 7V FFT = 128 (purple), 7V FFT = 256 (blue), N FFT = 350 
(red) and jV FFT = 375 (green). 



ures up to half the Nyquist frequency for our default choice 
of AWr = 350 3 , k nq /2 = ttAW /(2L box ) = 0.37/iMpc" 1 
which is beyond the range typically used in BAO fitting 
when assuming linear theory. 

The left panel in Fig. [3] shows the ratio of the power 
spectra, Pgg, P&e and Pgg measured at z — 0, to the power 
spectra measured at z = 5 scaled using the ratio of the 
square of the linear growth factor at z — 5 and z — for 
ACDM. It is clear from this plot that all P(k) evolve as ex- 
pected in linear theory on the largest scales. Note a linear 
scale is used on the a>axis in this case. In the right panel in 
Fig. [3] all the power spectra have been divided by the linear 
theory matter power spectrum measured from the simula- 
tion at z = 5, scaled using the ratio of the linear growth 
factor at z = 5 and z = 0. T his removes the sampling v ari- 
ance from the plotted ratio (|Baugh fc Efstathioulll994 ). In 
both panels, the error bars represent the scatter over eight 
simulations in ACDM averaging the power spectra after im- 
posing the distortions along the x, y or z axis in turn. From 
this figure we can see that the non-linear velocity divergence 
power spectra can be substantially different from the mat- 
ter power spectrum on very large scales k ~ 0.03/iMpc -1 . 
The linear perturbation theory assumption that the velocity 
divergence power spectra is the same as the matter P(k) is 
not valid even on these large scales. In the case of ACDM 
this difference is ~ 20% at k = O.lftMpc -1 . Note in the 
right panel in Fig. [3] the 10% difference in the ratio of the 
cross power spectrum to the matter power spectrum, on the 
largest scale considered, indicates that we have a biased es- 
timator of 6 which is low by approximately 10%. 

We find that the Pgg and Pgg measured directly from 
the simulation differ from the matter power spect rum by 
more then was reported bv lPercival fc White! (|2009t) . These 
authors did not measure Pgg and Pgg directly, but instead 
obtained these quantities by fitting Eq. [13] to the redshift 
space monopole power spectrum measured from the simu- 
lations. In Fig. [4] we plot the same ratios as shown in the 
right panel of Fig. [3] measured from one ACDM (left panel) 
and SUGRA (right panel) simulation. From our simulations 
it is possible to find a realisation of the density and velocity 
fields where the measured matter power spectrum and the 
velocity divergence power spectra are similar on large scales. 



the matter perturbation and r is the conformal time, dt = 
a(r)dr, then the linear continuity equation becomes 

9 = V -u = -aHfS . (12) 

Throughout this paper we normalise the velocity diver- 
gence as 6(k, a)/[— aH(a)f(£lin,(a), 7)], so 6 — 8 in the lin- 
ear regime. The volume weighted velocity divergence power 
spectrum is calculat ed from the simulati ons according to the 
prescription given in lScoccimarrol (|2004 ). We interpolate the 
velocites and the densities onto a grid of 350 3 points and 
then measure the ratio of the interpolated momentum to 
the interpolated density field. In this way, we avoid having 
to correct for the CIC assignment scheme. A larger grid di- 
mension could result in empty cells where S — >• 0. A FFT 
grid of 350 3 was used to ensure all grid points had non- 
zero density and hence a well defined velocity at each point. 
We only plot the velocity power spectra in each of the fig- 



Having found that the measured Pgg and Pgg differ 
significantly from Pgg, we now test if the grid assignment 
scheme has any impact on our results. As explained in Sec- 
tion 14.21 the velocity P(k) are computed by taking the 
Fourier transform of the momentum field divided by the 
density field to reduce the impact of the grid assignment 
scheme ()Scoccimarrd 120041 ) . IPueblas fc Scoccimarrol ()2009t ) 
showed that the CIC assignment scheme affects the mea- 
sured P(k) beyond ~20% of the Nyquist frequency. In Fig. 
[5]we show the power spectrum measurements for four dif- 
ferent FFT dimensions to show the scales at which we get 
a robust measurement. For Af FFT = 350 the power spectra 
have converged on scales up to k ~ 0.2/iMpc . 

In the top row of Fig. [6] the ratios Po(k)/P r (k) and 
Pi{k)/Po(jk) are plotted as solid lines in the left and right 
hand panels respectively. In this figure we have overplotted 
as grey dashed lines, the ratio of the redshift space monopole 
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Figure 6. The left hand column shows the ratio of the monopole of redshift power spectra to the real space power spectra at z = 
and 2 = 1. The right hand column shows the ratio of the quadrupole to monopole moment of the redshift space power spectra at z = 
and 2 = 1. Different rows show different dark energy models as labelled. Top row: The ratio of the redshift and real space power spectra 
in ACDM are plotted as solid lines in the left panel. The dashed lines represent the same ratio using Eq. 1131 for the monopole of the 
redshift space power spectrum. The dot-dash line represents the model given in Eg. I1UI which includes velocity dispersion effects. In the 
right panel the ratio of the quadrupole to monopole moment of the redshift space power spectra in ACDM are plotted as solid lines. The 
same ratio using Eg. 1141 for the redshift space power spectrum is plotted as dashed lines. Middle row: Same as the top row but for the 
SUGRA quintessence model. Bottom row: Same as the middle row but for the CNR quintessence model. 



moment to the real space power spectrum where 

Po(k) = Pss(k) + rJPso(k) + \fPeo{k) . (13) 
3 5 

On scales 0.05 < fe(^Mpc~ 1 ) < 0.2, this model for the 
redshift space power spectrum reproduces the measured 
P a (k, fi) and is a significant improvement compared to Eq. 
[U This form does not include any modelling of the damp- 
ing due t o velocity d is persio n. The extended model pro- 
posed by IScoccimarrol f20o3) given in Eq. [101 which does 
include damping, is also plotted as a black dot-dashed line 
for ACDM in the top row in Fig. [6] The redshift space 
quadrupole to monopole ratio in the quasi-linear regime, 



including the velocity divergence power spectra, is 

Ei = IE™ + SfPoo {m 
P S Pea + if Pee + \f 2 Pee ' 

This model does well at reproducing the ratio of the redshift 
space to real space power spectrum, although it underpre- 
dicts the ratio on scales k < 0.02ftMpc . The correspond- 
ing plots for the SUGRA and CNR models are shown in the 
middle and bottom rows of Fig. [6] It is clear that includ- 
ing the velocity divergence power spectrum in the model for 
Po an d Pi, produces a good fit to the measured redshift 
space power in both quintessence models on scales up to 
k ~ 0.2/iMpc" 1 . 



10 E. Jennings, C. M. Baugh, S. Pascoli 




P es (Mpc/h) 3 P, 5 (Mpc/h) 3 

Figure 7. Non-linear velocity divergence auto (left) and cross (right) power spectrum plotted as a function of the non-linear matter 
power spectrum at z = 0, 1 and 2 in three quintessence models and ACDM, as labelled. The ratio of the velocity divergence power spectra 
to the matter power spectrum at each redshift is plotted in the smaller panels beneath each main panel. 



5 RESULTS II: THE DENSITY VELOCITY 
RELATION 

In Section [5. II we examine the relationship between the non- 
linear matter and velocity divergence power spectra in differ- 
ent cosmologies. In Section [572] we study the redshift depen- 
dence of this relationship and provide a prescription which 
can be followed to generate predictions for the non-linear 
velocity divergence power spectrum at a given redshift. 

5.1 Dependence on cosmological model 

The linear continuity equation, Eq. 1121 gives a one to 
one correspondence between the velocity and density fields 
with a cosmology dependent factor, /(fi m ,7). Once the 
overd ensities become non -linear, this relationship no longer 
holds. iBernardeaul l| 19921 ) derived the non-linear relation be- 
tween S and 8 in the case of an initially Gaussian field. 
IChodorowski fc Lokasl (|l997T j extended this relation into the 
weakly non-linear regime up to third order in perturbation 
theory and found th e result to be a third order po lynomial 
in 6. More recently, iBilicki fc Chodorowskl (|2008l ) found a 
relation between 9 and 5 using the spherical collapse model. 
In all of these relations, the dependence o n cosmological pa- 
rameters was found t o be extremely weak (jBernardeaull 19921 : 
iBouchet et~al]|l995l ). The velocity divergence depends on 
O m and f^A, in a standar d ACDM cosmology, only through 
the linear growth rate, / l)Scoccimarro et aDll999l ). 

We showed in the previous section that including the 
velocity divergence auto and cross power spectrum accu- 
rately reproduces the redshift space power spectrum for a 
range of dark energy models on scales where the Kaiser for- 
mula fails. The quantities in Eqs. [14] and \W\ can be cal- 
culated if we exploit the relationship between the velocity 
and density field. In Fig. [7] we plot the velocity divergence 



auto (left panel) and cross (right panel) power spectrum as 
a function of the matter power spectrum for ACDM and 
the three quintessence dark energy models. We find that the 
density velocity relationship is very similar for each model 
at the redshifts considered, with only a slight difference for 
the SUGRA model at high redshifts and at small scales. The 
departure of the SUGRA model from the general density ve- 
locity relation is due to shot noise, which affects the power 
spectrum most at these scales in the SUGRA model as it 
has the lowest amplitude. We have verified that this effect is 
due to shot noise by sampling half the particles in the same 
volume, thereby doubling the shot noise, and repeating the 
P(k) measurement to find an even larger departure. Fig. [7] 
shows the independence of the density velocity relation not 
only of the valu es of cosmol o gical parameters, as found in 
previous works, IBernardeaul (|l992T ). but also a lack of de- 
pendence on the cosmological expansion history and initial 
power spectrum. 

Fitting over the range 0.01 < k(h/Mpc< 0.3), we 
find the following function accurately describes the rela- 
tion between the non-linear velocity divergence and mat- 
ter power spectrum at z = to better than 5% on scales 
k < 0.3/iMpc" 1 , 



Pxy(k) = g(P ss (k)) 



<Xo 

OL2 + a 3 P S s(k) 



(15) 



where Pss is the non-linear matter power spectrum. For 
the cross power spectrum P xy = P$g, Qo = —12288.7, 
Qi = 1.43, Q2 = 1367.7 and a 3 = 1.54 and for P xy = Pge, 
q = -12462.1, qi = 0.839, a 2 = 1446.6 and a 3 = 0.806; 
all points were weighted equally in the fit and the units for 
a>o, cti and Q3 are (Mpc//i) 3/2 , (Mpc//t)~ 3 and (Mpc/ft)" 3 
respectively. The power spectra used for this fit are the av- 
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Figure 8. A schematic illustration showing how the z = non- 
linear matter power spectrum can be rescaled to find the velocity 
divergence power spectrum at any redshift z = z' . The upper two 
curves represent the non-linear matter power spectrum, Pi, in 
grey and the velocity divergence power spectrum, P2, plotted as a 
blue dashed line, at z = 0. The power in the first bin is represented 
as a filled circle for each spectrum. The lower two curves, P[ and 
Pj, are the non-linear matter and velocity divergence spectra at 
z = z' . The power in the first bin is represented as a filled triangle 
in each case. The fitting formula for g(Pi) (Eg. I15|l generates the 
non-linear velocity divergence power spectra at 2 = 0. Using the 
function given in Eq. ll7l the matter power spectrum Pi and (/(Pi) 
can be rescaled to an earlier redshift. The power in the first bin 
from the rescaled Pi and g(P\) are shown as an empty grey and 
blue circle respectively. Note that Pi and P2 have been artifically 
separated for clarity. 



erage P e , 
tions. 



Pse and Pss measured from eight ACDM simula- 



5.2 Approximate formula for P<s, 
arbitrary redshift 



and Pag for 



In perturbation theory, the solution for the density con- 
trast is expanded as a se ries around the background value. 
IScoccimarro et ail (1 19981 ) found the following solutions for S 
and 8 to arbitrary order in perturbation theory, 



<5(fc,r) = ^D n (r)5 n (fc) 

n=l 
00 

6(k,T) = Y t En(T)e n (k), 



(16) 



where Si(k) and 9\(k) are linear in the initial density field, 
82 and 62 are quadrati c in the initial density field etc. 
IScoccimarro et ail (|l998l ) showed that using a simple ap- 
proximation to the equations of motion, /(fi m ) = Sim, the 
equations become separable and E n (r) — D n (r) = D(r) n , 
where D(t) is the linear growth factor of density perturba- 
tions. We shall use these solutions for S(k,r) and 0(k,r) to 
approximate the redshift dependence of the density velocity 
relation found in Section [5T] This relation does not depend 



on the cosmological model but we shall assume a ACDM 
cosmology and find the approximate redshift dependence as 
a function of the ACDM linear growth factor. 

The fitting function given in Eq. [15] generates the non- 
linear velocity divergence power spectrum, Pse or Pee from 
the non-linear matter power spectrum, Pss at z = 0. Fig. [5] 
shows a simple illustration of how the function g(Pss) and 
Pss at z = can be rescaled to give the velocity divergence 
power spectra at a higher redshift, z' . Using the simplified 
notation in the diagram, where Pi = Pss, and given the 
function g(Pss), we can find a redshift dependent function, 
c(z), with which to rescale g(Pss(z = 0)) to the velocity di- 
vergence P(k) at z' . At the higher redshift, z' , the non-linear 
matter and velocity divergence power spectra are denoted as 
P[ and P'2 respectively in Fig. [H] 

Using the solutions in Eg. 1161 to third order in perturba- 
tion theory, see Appendix fX] we assume a simple expansion 
with respect to the initial density field, to find the following 



• Pi(z 

0,2') 



z ) which can 
> 5 (Pi)/c 2 ( Z = 



ansatz for the mapping P{(z — z') - 
be approximated as Pi (z = 0)/c 2 (z 
0, 2') where 

D(z) + D 2 (z) + D 3 (z) 
{ ' ' D(z') + D 2 (z') + D^(z') ' 

and D(z) is the linear growth factor. The equivalence of 
these mappings gives P1 — P2 = (Pi — g(Pi))/c 2 which allows 
us to calculate P2 at z = z' if we have Pi(z = 0), g(Pi(z = 0) 
and P[(z — z'). Writing this now in terms of Pss, instead of 
Pi, we have the following equation 

p fu 's _ g(Pss(k,z = 0))-P ss (k,z = 0) 



(17) 



c 2 (z 



0, 



+Pss(k,z') 



(18) 



where g(Pss) is the function in Eg. I15l and P xy is either the 
nonlinear cross or auto power spectrum, Pse or Pss. 

In the left panel of Fig. [5] we plot the ACDM non-linear 
power spectrum Pee at z = 0, 1, 2 and 3. The function given 
in Eq.[l8]is also plotted as red dashed lines using the factor 
c(z, z') given in Eq. [T7]and the ACDM linear growth factor 
at redshift z = 0, 1, 2 and 3 respectively. The ratio plot shows 
the difference between the exact Pee power spectrum and the 
function given in Eq. 1181 The right panel in Fig. [5] shows a 
similar plot for the Pse power spectrum. In both cases we 
find very good agreement between the scaled fitting formula 
and the measured power spectrum. Scaling the z = power 
spectra using this approximation in Eq. [T7] reproduces the 
non-linear z = 1,2 and 3, Pse to ~ 5% and Pee to better 
than 5% on scales 0.05 < fe(/iMpc _1 ) < 0.2. It is remarkable 
that scaling the z = fitting formula using c in Eg. 1 171 works 
so well at the different redshifts up to k < 0.3/i/Mpc and is 
completely independent of scale. 

To summarise the results of this section we have found 
that the quadrupole to monopole ratio given in Eq. [14] and 
the model in Eq. 1101 which includes the non-linear matter 
and velocity divergence power spectra at a given redshift 
z', can be simplified by using the following prescription. As- 
suming a cosmology with a given linear theory matter power 
spectrum we can compute the non-linear matter P(k) at 
z=0 and at the required redshift, z' , using, for example, the 
phenomenological mode l HALQFIT dSmith et al.l [2003h or 
the method proposed bv lCasarini et al" I (|2009h in the case of 
quintessence dark energy. These power spectra can then be 
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Figure 9. Non-linear velocity divergence auto and cross power spectrum, in the left and right panels respectively, measured from the 
ACDM simulations at z = (open grey squares), z = 1 (purple crosses), z = 2 (blue stars) and z = 3 (cyan diamonds). Overplotted as 
red dashed lines is the function given in Eg. I18l at redshifts z = 1,2 and 3. The lower panels show the function in Eq. [18] divided by the 
measured spectra at z = 1, 2 and 3. 



used in Eq. [18] together with the function g, given in Eg. 1151 
and the linear theory growth factor between redshift z = 
and z = z' to find the velocity divergence auto or cross power 
spectrum. As can be seen from Fig. [9] the function given in 
Eq. [18] agrees with the measured non-linear velocity diver- 
gence power spectrum to ~ 10% for k < 0.3/tMpc -1 and to 
< 5% for k < 0.2/iMpc" 1 for ACDM. We have verifed that 
this prescription also reproduces Pse and Pgg to an accuracy 
of 10% for k < 0.3/iMpc" 1 for the CNR, SUGRA and 2EXP 
models using the corresponding matter power spectrum and 
linear growth factor for each model. This procedure simpli- 
fies the redshift space power spectrum in Eq. [10] and the 
quadrupole to monopole ratio given in Eq. 1141 For the dark 
energy models considered in this paper, this ratio provides 
an improved fit to the redshift space P(k, fi) compared to 
the Kaiser formula and incorporating the density velocity 
relation eliminates any new parameters which need to be 
measured separately and may depend on the cosmological 
model. 



6 CONCLUSIONS AND SUMMARY 

One of the primary goals of future galaxy redshift surveys is 
to determine the physics behind the accelerating expansion 
of the Universe by making an accurate measurement of th e 
growth rate, /, of large scale structure (|Cimatti et al.l f2009). 
Measuring the growth rate with an error of less than 10% 
is one of the main science goals of Euclid, as this will allow 
us to distinguish modified gravity from dark energy models. 
With an independent measurement of the expansion history, 
the predicted growth rate for a dark energy model would 
agree with the observed value of / if general relativity holds. 



We use simulations of three quintessence dark en- 
ergy models which have different expansion histories, lin- 
ear growth rates and power spect r a com pared to ACDM. In 
a previous paper. [Jennings et al.l |2010f), we carried out the 
first fully consistent N-body simulations of quintessence dark 
energy, taking into account different expansion histories, lin- 
ear theory power spectra and best fitting cosmological pa- 
rameters fl m , fib and Ho, for each model. In this paper we 
examine the redshift space distortions in the SUGRA, CNR 
and 2EXP quintessence models. These models are represen- 
tative of a broader class of quintessence models which have 
different growth histories and dark energy densities at early 
times compared to ACDM. In particular the SUGRA model 
has a linear growth rate that differs from ACDM by ~ 20% 
at z = 5 and the CNR model has high levels of dark energy 
at early times, S1 D e ~ 0.03 at 2 ~ 200. The 2EXP model 
has a similar expansion history to ACDM at low redshifts, 
z < 5, despite having a dynamical equation of state for the 
dark energy comp onent. For more details on each of the dark 
energy models see I Jennings et all (|2010T ) . 

Redshift space distortions observed in galaxy surveys 
are the result of peculiar velocities which are coherent on 
large scales, leading to a boost in the observed redshift space 
power spectru m compared to the real space power spectrum 
|Kaiser]|l987l ). On small scales these peculiar velocities are 
incoherent and give rise to a damping in the ratio of the 
redshift to real space power spectrum. The Kaiser formula 
is a prediction of the boost in this ratio on very large scales, 
where the growth is assumed to be linear, and can be ex- 
pressed as a function of the linear growth rate and bias, 
neglecting all non-linear contributions. 

In previous work, using N-body si mulations in a pe - 
riodic cube of 300/i _1 Mpc on a side, ICole et all (|l994T ) 
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found that the measured value of = f/b, where b is 
the linear bias, deviates from the Kaiser formula on wave- 
lengths of 50/T 1 Mpc or more as a result of these non- 
linearities. lHatton fc Cold (1 19981 ) extended this analysis to 
slightly larger scales using the Zel'dovich approximation 
combined with a dispersion model where non-linear veloci- 
ties are treated as random perturbations to the linear the- 
ory velocity. These previous studies do not provide an accu- 
rate description of the non-linearities in the velocity field for 
two reasons. Firstly, the Zel'dovich approximation does not 
model the velocities correctly, as it only treats part of the 
bulk motions. Secondly, in a computational box of length 
300fo -1 Mpc, the power which determines the bulk flows has 
not converged. In this work we use a large computational 
box of side 1500ft _1 Mpc, which allows us to measure red- 
shift space distortions on large scales to far greater accuracy 
than in previous work. 

In this paper we find that the ratio of the monopole of 
the redshift space power spectrum to the real space power 
spectrum agrees with the linear theory Kaiser formula only 
on extremely large scales k < 0.03/iMpc -1 in both ACDM 
and the quintessence dark energy models. We still find sig- 
nificant scatter between choosing different axes as the line 
of sight, even though we have used a much larger simula- 
tion box than that employed in previous studies. As a result 
we average over the three power spectra, assuming the dis- 
tortions lie along the x, y and z directions in turn, for the 
redshift space power spectrum in this paper. Instead of using 
the measured matter power spec trum in rea l space , we find 
that the estimator suggested by I Cole et alj (| 19941 ). involv- 
ing the ratio of the quadrupole to monopole redshift space 
power spectrum, works better than using the monopole and 
agrees with the expected linear theory on slightly smaller 
scales k < 0.07/iMpc" 1 at z = for both ACDM and the 
quintessence models. 

As the measured redshift space distortions only agree 
with the Kaiser formula on scales k < 0.07/iMpc -1 , it 
is clear that the linear approximation is not correct on 
scales which are normally considered to be in the 'linear 
regime', k < 0.2MVIpc -1 . In linear theory, the velocity di- 
vergence power spectrum is simply a product of the mat- 
ter power spectrum and the square of the linear growth 
rate. In this work we have demonstrated that non- linear 
terms in the velocity divergence power spectrum persist on 
scales 0.04 < fcjTiMpc -1 ) < 0.2. These results agree with 
IScoccimarrd 1I2OO4Y who also found significant non-linear 
corrections due to the evolution of the velocity fields on large 
scales, assuming a ACDM cosmology. We have shown that 
including the non-linear velocity divergence auto and cross 
power spectrum in the expression for the redshift space P(k) 
leads to a significant improvement when trying to match the 
measured quadrupole to monopole ratio for both ACDM and 
quintessence dark energy models. 

Including the non-linear velocity divergence cross and 
auto power spectra in the expression for the redshift space 
power spectrum increases the number of parameters needed 
and depends on the cosmological model that is used. Using 
the non-linear matter and velocity divergence power spectra 
we have found a density velocity relation which is model 
independent over a range of redshifts. Using this relation it 
is possible to write the non-linear velocity divergence auto 
or cross power spectrum at a given redshift, z' , in terms 



of the non-linear matter power spectrum and linear growth 
factor at z = and z = z . This formula is given in Eq.ll8lin 
Section l5.2l We find that this formula accurately reproduces 
the non-linear velocity divergence P(k) to within 10% for 
k < 0.3/iMpc" 1 and to better than 5% for k < 0.2/iMpc" 1 
for both ACDM and the dark energy models used in this 
paper. 

It is clear that including the non-linear velocity diver- 
gence terms results in an improved model for redshift space 
distortions on scales k < 0.2/iMpc" 1 for different cosmologi- 
cal models. Current galaxy redshift su rveys can provide onl y 
very weak constraints on P se and Pee (jTegmark et al.ll2002h . 
The relation given in this paper between the non-linear ve- 
locity divergence and matter power spectra will be useful for 
analysing redshift space distortions in future galaxy surveys 
as it removes the need to use noiser and sparser velocity 
data. 
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APPENDIX A: APPROXIMATE FORMULA 
FOR P se AND P ge FOR ARBITRARY REDSHIFT 

Eq. 18 in this paper relates P xy (z') — Pss(z') at z = z' to 
the same expression at redshift z = using a variable c 2 . 
Note from Eq. 15 g(Pss(z = 0)) = P xy {z = 0) in Eq. 18. 
From Eqs. 16 in our paper and using the result by Scoc- 
cimarro et al. 1998 we can write the following solutions 
for 9 and 8 in terms of scalings of the initial density field 
l|Bernardeau et alj|2002h . 

e(z) = D{ z )e 1 + D 2 (z)e 2 + D 3 (z)9 s + --- (Ai) 



and 

5{z) = D(z)S 1 + D 2 (z)8 2 + D 3 {z)S 3 + - 



(A2) 



Squaring these expressions and emsemble averaging we can 
write the velocity divergence power spectrum and the matter 
power spectrum to third order in perturbation theory as 

Pgg(z') ~ <\D{z')6 1 + D 2 (z')8 2 + D 3 (z')8 3 \ 2 > (A3) 
Pss(z') ~ <\D(z')5 1 +D 2 (z')S 2 + D 3 (z')S :i \ 2 > . (AA) 

Using the fact that \D0i + D 2 9 2 + D 3 6 3 \ < \D9i \ + \D 2 6 2 \ + 
\D 3 9 3 \ we can approximate this as 

Peo(z') « {Diz')^] + D 2 {z')\9 2 \ + D 3 {z')\9 B \) 2 > (A5) 
Pss(z') ^< (D(z')\5i\ + D 2 (z')\S 2 \ + D 3 (z')\S 3 \) 2 > , (A6) 
and we assume that 

< \D(z')6 1 + D 2 (z')6 2 + D 3 (z')9 3 \ 2 > (A7) 

- < (£>(z')|0i| +D 2 (z')\e 2 \ +D 3 (z')\e 3 \) 2 >~ 

< \D(z')Si + D 2 (z')S 2 + D 3 (z')S 3 \ 2 > 

- < {D(J)\6i\ +D 2 (z')\S 2 \ + D 3 (z')\S 3 \) 2 > . 
Taking the difference of the two power spectra we have 

Poe(z') -Pss(z') ~ 

< (£>(*')l*i| + D 2 {z')\6 2 \ +D 3 (z')\e 3 \) 2 > 

- < (D(z')\di\ + D 2 (z')\5 2 \ + D 3 {z')\8 3 \) 2 > (A8) 

and as x 2 — y 2 = (x — y) (x + y) we can rewrite this as 

Pee{z) -Pss{z') - (A9) 

< [DQOxl - \Si\) + D 2 (\9 2 \ - \5 2 \) + D 3 (\9 3 \ - \5 3 \)} 
x[D(\e 1 \ + \8 1 \) + D 2 (\9 2 \ + \8 2 \) + D 3 (\e 3 \ + \8 3 \)] > . 

Multiplying out the rhs of this equation and denoting the 
modulus of variable |a;| as i for simplicity, we have 

Peo(z') - Pss(z') - 

< {D 2 [6l - 5l] + D 3 [(0i -5^(62 + 82) + {9x +8^(62-8 

+D i [{e 1 - Si){e 3 + s 3 ) + (el - s 2 2 ) + (Ox + 5i)(6> 3 - s 3 )] 

+D S [{9 2 - S 2 )(6 3 + 8 3 ) + (62 + S 2 ){9 3 - &)] + D a [6 2 3 - <5| 
and then taking out a factor of [8\ — 8 2 ] on the rhs we have 
P ee (z) - P ss (z') ~ (All) 



< [0\ — 8 2 }{D + D 



,62 + 82 , 62-82 



+D*[ 



+D 5 \2 



; r <->s + 8 3 9 2 — 8 2 



+ 



1+S1 di- 81 
3 — 8 3 , 



t + h ei-8* e 1 -s 1 

" "2 - ^3*2, , , v " ; -», 



- Ss8 2 ] fir Oj - 8 2 

91-81 J + [ e 2 -s'( 



]}> 
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As 6*i and Si are linear in the initial density contrast, which 
we assume to be different to the linear density contrast, 
9i ~ A ~ Si and 62 ~ £2 ~ Si + S 2 is quadratic in the 
initial density contrast and #3 ~ 83 ~ Si + 5 2 + 8f is cubic 
in the initial density field, we assume 61+62 ~ 5i + 82, 
61 + 6.i ~ 5i + £3 and #1 — #2 ~ <5i — #2, #1 — 6*3 ~ Si — S3 so 
the fractions in the above equation are unity and 

Peo(z') - Pst(z') (A12) 
- < [6<i — Si] > {D 2 + 2D [i + 3D 4 + 2D 5 + D G } 
<[6'i- 5!] > {D(z') + D 2 (z') + D a (z')} 2 

Similarly for Pee (z) — Pss (z) we have 

Pee(z) - Pss(z) (A13) 
~ <[6 2 -S 2 ]>{D(z) + D 2 (z) + D 3 (z)} 2 

Taking the ratio of the two previous equations, the redshift 
independent factor [6 2 — S 2 ] cancels and we obtain the fol- 
lowing ansatz 



[D(z') + D(z') 2 + D{z'ff 
[D(z) + D(z) 2 + D(z) 3 ] 2 



(A14) 



Pee(z') - Pss(z') 
Pee{z) - Pss(z) 

which is the expression in Eq. 18 in the paper for z = 0. A 
similar approximation works for the cross power spectrum 
Pie- 
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